Quantum Circuit Design for Solving Linear Systems of Equations 
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Recently, it has been demonstrated that quantum computers can be used for solving linear systems 
of algebraic equations with exponential speedup compared with classical computers. Here, we 
present an efficient and generic quantum circuit design for implementing the algorithm for solving 
linear systems. In particular, we show the detailed construction of a quantum circuit which solves a 
4x4 linear system with 7 qubits. It consists of only the basic quantum gates that can be realized with 
present physical devices, implying great possibility for experimental implementation. Furthermore, 
the performance of the circuit is numerically simulated and its ability to solve the intended linear 
system is verified. 



Quantum computers are devices that take di- 
rect advantage of quantum mechanical phenom- 
ena such as superposition and entanglement to 
perform computations [1]. Because they com- 
pute in ways that classical computers cannot, 
for certain problems quantum algorithms pro- 
vide exponential speedups over their classical 
counterparts. For example, in solving problems 
related to factoring large numbers [2] and simu- 
lation of quantum systems [51-[T7] , quantum al- 
gorithms are able to find the answer exponen- 
tially faster than classical algorithms. Recently, 
A. W. Harrow, A. Hassidim and S. Lloyd dHj 
proposed a quantum algorithm for solving linear 
systems of equations with exponential speedup 
over the best known classical algorithms. 

With the theoretical potential of the algo- 
rithm, an equally important issue is to ren- 
der the algorithm viable for experimental imple- 
mentation. Therefore, in this work we present 
a generic quantum circuit design based on the 
theoretical framework outlined by the algo- 
rithm. Particularly, the original work [TS] did 
not elaborate on the circuit design for imple- 
menting A~^-controlled rotation when only Xj is 
available. Hence in this work we present a scal- 
able design of the quantum circuit for finding 
the reciprocals of the eigenvalues Xj and stored 
them in a quantum register with a superposi- 
tion of states that encode the values of XJ^ in 



binary format. The runtime scaling of this in- 
verting circuit is O(logA^), which retains the ad- 
vantageous O(logA^) runtime for the overall al- 
gorithm. For the purpose of showing a concrete 
numerical example of how this scheme works, 
we choose a specific linear system (described by 
a matrix A of dimension 4 x 4) as the numerical 
example for of the algorithm. Since our circuit 
involves only 7 qubits and is composed of only 
basic quantum gates, our work bridges the the- 
oretical development of the algorithm with the 
possibility of physical implementation by exper- 
imentalists. 

The algorithm [TS] solves the problem Ax = b 
where A, a Hermitian d-sparse NxN matrix 
and b, a unit vector are given. The major steps 
of the algorithm can be summarized as the fol- 
lowing: (1) Represent the vector 6 as a quantum 
state \b) = X]i3=i^jN) stored in a quantum reg- 
ister (termed B), where \i) represents a basis 
state of the register B where the states of the 
qubits represent the binary coding of i. In a sep- 
arate quantum register (termed C) of t qubits, 
initialize the qubits by transforming the regis- 
ter to state |r) from |0). Here |r) repre- 
sents a basis state of the register C where the 
states of the qubits represent the binary cod- 
ing of r; (2) Apply the conditional Hamiltonian 
evolution \T){Tf (E> e'-^rto/T. (3) ^pp^y 

the quantum inverse Fourier transform to the 
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FIG. 1: Overall circuit design. The registers 

from top to bottom are respectively the 
ancillary qubit, register L, M, C and B. The 
initial state is denoted as 1 00006). 
IV;/) = ^1- f^|0) + W represents 

uncomputation. 



register C. Denote the basis states after quan- 
tum Fourier transform as |fc). At this stage, 
the amplitudes of the basis states are concen- 
trated on k values that satisfy ~ where 
Afe is the fc-th eigenvalue of the matrix A; (4) 
Add an ancilla qubit and apply conditional ro- 
tation on it, controlled by the register C with 
1^) ~ \^k)- The rotation tranforms the qubit to 
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^ 1 0) + r" 1 1) • This is a key step of the al- 
gorithm and it involves finding the reciprocal of 
the eigenvalue Xj quantum mechanically, which 
is not a trivial task on its own. We will elabo- 
rate on how to achieve this later in this work; (5) 
Uncompute the registers b and C, and (6) Mea- 
sure the ancilla bit. If it returns 1, the register 
b of the system is in the state X]j'=i Pj-^j^^l'^^j) 
up to a normalization factor, which is equal to 
the solution \x) of the linear system Ax = b. 
Here \uj) represents the j-th eigenvector of the 
matrix A and let \b) = J2^=i l^jl^j)- 

Based on the theoretical framework of the 
algorithm outlined above, here we present a 
generic quantum circuit design that implements 
the algorithm to find the solution \x) to the 
linear equation A\x) = \b). The circuit con- 
tains registers B, C, M, L and an ancillary 
qubit (Figure [I]). The detailed matrix forms of 
the quantum gates are presented in Appendix. 



Here the time parameter to for all Hamiltonian 
simulations is defined to be 2tt. Register B is 
used to store the value of \b). Register C, M 
and L contain t, m and I qubits respectively, 
the Walsh-Hadamard transform on Register C, 
the controlled Hamiltonian simulation and the 
inverse Fourier transform ideally gives a state 
T,j l^j\>'j)\uj) [H in Register B and C. At 
the same time, Walsh-Hadamard transform and 
controlled global phase shift R^z is applied to 
register L and M (Figure[T]) and the state of the 
two registers becomes X]p Z]s^^P(*^*o)|s)b), 
where \p) and \s) represent the basis states 
of the register M and L respectively. Details 
of this part of the circuit is shown in Figure 
[5l After the inverse Fourier transform is ex- 
ecuted on register C, we use its \Xj) states 
stored in register C as a control for a Hamil- 
tonian simulation exp(iifoto) that is applied on 
register Af (Details in Figure [6]) . Here Hq is 
defined as diagonal matrix whose diagonal ele- 
ments are (1, 2, 2™^^). Due to the unitary 
nature of exp(ji/oio), the operation exp(iifoio) 
can be readily decomposed into a quantum cir- 
cuit [m [20] that consists of only basic quantum 
gates. 

We further establish the control relationship 
between Register L and the exp(ii/oto) by us- 
ing the the fej-th qubit (Figure [6]) of the register 
L to control a exp[— ip(^ ^j-^io)] simulation 
that acts on register M. The values of binary 
numbers stored in register L are then able to 
determine the time parameter t in the over- 
all Hamiltonian simulation exp(— ii/oi). Fol- 
lowing such construction, after the Hamiltonian 
simulation exp(— iHgio) the state of the sys- 
tem is |0) (8) J2]=i T,l=o^ ELV ^jexp[i 2^io 
(2'-Ajs)] \s)\p) \\j)\uj). With to = 2n, the 
states Is) is concentrated on values of s = |- in 
the above expression (Denote such |s) states as 
||-)). This is because of the fact that sums 

of the form ^^^QGxp{2Trikj^) vanish unless 
r — OmodiV. Denote the components of the 
superposition states in register L as ||-). Ro- 
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tate the anciUa bit with the angle shift con- 
trolled by the ||-) states stored in register L: 



Ryii-m « TjlO) + Ijll) ^ith C be- 

ing a constant. Uncomputing the three reg- 
isters (Figure [T]), we are left with a state in 
register B with the ancilla qubit proportional 

to Ej=i [^l-C2/A2|0) + C/A,|l)] ^ |000) 

Measuring the ancilla qubit and if we ob- 
tain |1), the state of the system will collapse 
to 1 1000) '^'^^^i PjXj^luj), where \uj), where 
the first register represents the solution. If the 
ancilla bit measures |0), the system is back to 
its original state 1 00006). If the ancilla bit mea- 
sures |1), the state stored in the register B is 
proportional to the solution of the linear sys- 
tem of equations A\x) — \b) and the register 
B is no longer entangled with any other regis- 
ters in the system. Therefore we have obtained 
the solution of the linear system of equations in 
form of a quantum superposition state stored in 
the register B. 

There are mainly two factors that influence 
the number of qubits needed in register L and 
M. The first one is the desired bit precision for 
encoding both Xj and . Clearly if one desires 
to encode a number with precision e — 0(2''), 
at least q qubits are needed for the register, 
which is true for any generic binary registers. In 
addition to this factor, for the particular roles 
that the register L and M play in the circuit, 
the condition number k of A and the condition 
number of the matrix A is also important. Be- 
cause in order for the inversion circuit to work, 
2™ must be larger than any eigenvalue of A i.e. 
"^ > [log Xmax] ■ On the other hand in order 
to encode all values of t^- register L must have 
enough qubits such that even the largest one 
can be encoded i.e. I > [log Xmin^ ■ There- 
fore m — I log K is a general constraint that 
m and I must comply with. 

A detailed runtime analysis in [18; shows that 
the runtime scales as 0{\og{N)(P) where d is 
the sparsity of A. The circuit that we present 
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FIG. 2: The final state of the algorithm for 
solving the 4x4 system. The horizontal axis 
represents the decimal value that corresponds 
to a basis state of the seven-qubit system. The 
vertical axis represents the real parts of the 
probability amplitudes that corresponds to a 
certain basis. 



here adds to the algorithm the additional com- 
plexity of controlled Hamiltonian simulation 
exp(ii/oio)- However, it is shown in |21j that 
e^^* can be simulated in time 0{\og{N)(Pt). 
Therefore the Hamiltonian simulations in this 
circuit design only add at most a constant over- 
head to the cost scaling of the algorithm. 

With the generic circuit model outlined 
above, now we present a 7-qubit circuit for solv- 
ing a system with A of dimension 4 x 4 as a 
numerical example. For the numerical example 
we specifically choose the form of matrix A such 
that it allows for a very simple ad hoc design of 
the circuit for finding the reciprocal modulo 2* 
(in the 7-qubit case 2* = 16). For this example 
we choose 



A = 



/15 9 5 -3^ 

9 15 3 -5 

5 3 15 -9 

V-3 -5 -9 15/ 



(1) 



Hence A is a Hermitian matrix with the eigen- 
values Xi = 2*~^ and corresponding eigenvec- 
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FIG. 3: Simulation results on the dependence 
of final state on the value of the parameter r. 
The first barplot shows the probability 
distributions over the four basis states of 
Register M with different values of r in the 

controlled rotation gates acted upon the 
ancilla bit. The horizontal dashed lines show 
the analytical values of the probabilities. The 
two plots following show the probability of 
getting the solution as well as the fidelity 
{x'\x) as functions of r. 



tor \u,) = where rep- 

resents the state of register C which encode the 
number j in binary form, Sij is the Kronecker 
delta, and the index i runs from 1 to 4. Fur- 
thermore, we let b = 1,1,1)^. Therefore, 

1 6) = J2'j=il^jWj) ^^'^ ea.ch f3j = \. To com- 
pute the reciprocals of the eigenvalues, a quan- 
tum swap gates is used (Figure |4]) to exchange 
the values of the first and third qubit. Therefore 
there is no need for specific auxiliary registers 
M and L (Figure[T]) in this example. The eigen- 
values are inverted to their reciprocals with the 
swap gate in the following fashion. For exam- 
ple, the eigenvalue A4 =8 is encoded as |1000) 
in register C, after applying the swap gates 
it is transformed to |0010), which is equal to 
i X 2* = 2. 

o 

Figure |4] shows the circuit. We use Group 
Leader Optimization Algorithm to find the cir- 
cuit decomposition of the Hamiltonian simula- 



sition of e'^ 



tion operator e'"^ 1^ (see [111 [ID] for details, also 
see [U 1^ for general methods). Here we let 
to = 27r. The resulting quantum circuit for 
is shown in Figure |4| With the decompo- 
readily available, the operators 
e"'T, e°"'T and e^^^ can be implemented by 
simply multiplying the angle shifts in all the ro- 
tation gates in the circuit for e'"* 1^ by a factor 
of 2, 4, and 8 respectively [19]. 

The final state of system, conditioned 
on obtaining |1) in the ancilla bit, is 
8|mi) -|-4|?i2) +2|m3) + |U4))«)|0000)®|1). 

in the standard bases, it becomes 
(-|00) + 7|01) + 11|10) + 13|11)), which is 
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1 

/340 



proportional to the exact solution of the system 
X = ^(-1,7,11,13)^ (Figure B. The quan- 
tum circuits (Figure |4| outlinea above is then 
simulated and the statistics of final states that 
are computed from the simulations (Figure [3]). 

The exponent r in the rotation gates R{i^^) 
(Figure |4]) is an important parameter that de- 
termines the accuracy of the final state of the 
system, which contains the solution to the lin- 
ear system of equations. As can be seen in Fig- 
ure [3j when the value of r is sufficiently large, 
say larger than 6, the probabilities of measur- 
ing each basis state |00), |01), |10) and |11) con- 
verge to the analytical values. Those analytical 
values correspond to the analytical solution of 
the 4 X 4 system x = 7^ 13)^^ p^j._ 

thermore, the fidelity of the results {x'\x) in- 
creases as r is increased. However, in Figure [3j 
as r grows beyond around 5, the probability of 
measuring the ancilla bit as |1) decays, which 
indicates that as r is increased, the solutions 
obtained in the final state in register b becomes 
more accurate yet less probable to attain. 

In conclusion, in this work we present con- 
structions of quantum circuits, both general and 
particular, for implementing the algorithm for 
solving linear systems of equations. Our results 
may motivate experimentalists with the capa- 
bility of addressing 7 or more qubits and ex- 
ecute basic quantum gates on their setups to 
implement the algorithm and verify its results. 
Future work regarding the improvement upon 
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FIG. 4: Quantum Circuit for solving Ax — 6 with v44x4 being the matrix shown in equation 
The first qubit at the top of the figure is the ancilla bit. The four qubits in the middle stand for 
the register C. The two qubits at the bottom is the register that stores the vector h. The detailed 
matrix forms of the quantum gates are presented in Appendix. 



the algorithm needs to be focused on topics such 
as extending the algorithm to cases where the 
matrix A has a high condition number |23| . 
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APPENDIX 

The matrix representations of the quantum 
gates and algorithms used in this work are the 
following: X, Y and Z gates which are the Pauli 
operators a.j;, Cy and and Hadamard gate: 



S gate and T gate (or |) gate are 




Square root of X gate (or NOT gate) and its 
conjugate are 



1 + z 1 l + 

2\l-i l + ij' 2\l + i l-ij 

(5) 

Rotation gates are 
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FIG. 5: Detailed quantum circuit of the generic circuit implementation. This figure shows the 
Walsh-Hadamard transform, |r) -controlled Hamiltonian simulation e*"^* and Fourier transform on 

the register C. 
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